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Abstract 

We present a review of our previous numerical studies on non equilibrium 
vortex dynamics in Josephson Junction arrays (JJA) driven by a dc current. 
Dynamical phase diagrams for different magnetic fields, current directions 
and varying temperature are discussed and compared. 

First, the effect of thermal fluctuations in a current driven diluted vortex 
lattice (VL) is analyzed. The case of / = 1/25 (where / is the fraction of flux 
quanta per plaquette in the array) is considered and the phase diagram as a 
function of the driving current and temperature is analyzed. In equilibrium, 
this system has a weakly first-order melting transition of the vortex lattice, 
which coincides with a depinning transition. When a low current is applied, 
the "longitudinal" depinning transition occurs at a temperature lower than 
the melting transition. More interestingly, for large currents (well above the 
critical current) there is an analogous sequence of transitions but for the 
transverse response of a fast moving VL. There is a transverse depinning 
temperature below the melting transition of the moving VL. 
We also discuss the dependence with the direction of the applied dc current 
of the transport properties of diluted vortex arrays on a square JJA at low 
temperatures. We show that orientational pinning phenomenon leads to 
a finite transverse critical current when the bias current is applied in the 
directions of high symmetry and it leads to an anomalous transverse voltage 
when vortices are driven away from the favorable directions. In addition, 
the effect of disorder in the transport properties of square JJA with a dc 
current applied in the "diagonal direction" ([11] direction) is analyzed and 
a finite transverse voltage is also observed in this case. 
The case of a fully frustrated square JJA, corresponding to / = 1/2, driven 
by a dc current and with thermal fluctuations is also discussed. In equilib- 
rium, the low temperature phase has two broken symmetries: the U(l) sym- 
metry, corresponding to superconducting coherence, and the Z2 symmetry 
corresponding to the periodic order of the VL, which forms a "checkerboard 
pattern". At high currents (well above the critical current) two well sepa- 
rated transitions are observed. The order of the checkerboard vortex lattice 
(discrete Z2 symmetry) is destroyed at a much lower temperature than the 
transverse superconducting coherence (continuous U(l) symmetry). 

Keywords: Josephson, vortex, phase transitions, non-equilibrium. 

1.1 INTRODUCTION 

The behavior of superconducting vortices in the presence of periodic pinning 
shows very rich static and dynamic phenomena ^0311111110111111013113 E3 
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tion between the repulsive vortex-vortex interaction and the attractive peri- 
odic pinning potential results in novel vortex structures at low temperatures 
[SI El El Ell El EEZ] • The equilibrium phase transitions of these vortex struc- 
tures and their various dynamical regimes when driven out of equilibrium are 
of great interest both experimentally [H I^J SI |31 El El El El H21 and theoreti- 
cally [SliainiCEaCnaElIiaCEia^ Sev- 
eral techniques have been developed to fabricate in superconducting samples 
an artificial periodic pinning structure: thickness modulated superconduct- 
ing films [J, superconducting wire networks |2j, Josephson junction arrays 
[3 IU [3] , magnetic dot arrays [H] , submicron hole lattices El , and pinning 
induced by Bitter decoration [9j- The ground states of these systems can 
be either commensurate or incommensurate vortex structures depending on 
the vortex density (i.e the magnetic field). In the commensurate case, a 
"matching" field is defined when the number of vortices N v is an integer 
multiple of the number of pinning sites N p : N v = nN p . A "submatching" or 
"fractional" field is defined when N v is a rational multiple of N p : N v = fN p 
with / = p/q. One of the main properties of periodic pinning is that there 
are enhanced critical currents and resistance minima both for fractional and 
for matching magnetic fields, for which the vortex lattice is strongly pinned. 

In the case of JJA |3] the discrete lattice structure of Josephson junctions 
induces an effective periodic pinning potential (the so-called "egg-carton" 
potential) which at low temperatures confines the vortices at the centers 
of the unit cells of the network [I]. There are strong commensurability 
effects for submatching fields / = p/q, for which the vortices arrange in an 
ordered qxq superlattice that is commensurate with the underlying array of 
junctions. The transition temperature T c (f) and the critical current I c (f) 
have maxima for rational / = p/q, which have been observed experimentally 
[SI El- Moreover, as we will discuss here, the model that describes the physics 
of the JJA can be thought as a discrete lattice London model for thin film 
superconductors with periodic arrays of holes. However, this comparison 
can be valid only for low submatching fields since it can not describe the 
effects of interstitial vortices. 

The equilibrium phase transitions at finite temperatures of two dimen- 
sional systems with periodic pinning have been studied in the past jl H 1121 
ESElEil- It is possible to have a depinning phase transition of the com- 
mensurate ground states at a temperature T p and a melting transition of the 
vortex lattice at a temperature Tm jllM12lH^lll41H5j . Franz and Teitel |14j 
have studied this problem for the case of submatching fields. For < T < T p 
there is a pinned phase in which the vortex lattice (VL) is pinned commensu- 
rably to the periodic potential and has long-range order. For T p < T < Tm 
there is a floating VL which is depinned and has quasi-long-range order. 
For high submatching fields (/ > 1/30 ) both transitions coincide, T p = Tm, 
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while for low submatching fields (/ < 1/30) both transitions are different 
with T p <T M HaHSl. 

The non-equilibrium dynamics of driven vortex lattices interacting either 
with random or periodic pinning shows an interesting variety of behavior 
d 1201 HH I22H231 HI HSHm US EBB ED El E31E1ES]. Many recent studies 
have concentrated in the problem of the driven VL in the presence of random 
pinning [213 E01 EU E21 EH1 EH EH] • When there is a large driving current the 
effect of the pinning potential is reduced, and the nature of the fastly moving 
vortex structure has been under active discussion. The moving vortex phase 
has been proposed to be either a crystalline structure, a moving glass, a 
moving smectic or a moving transverse glass [291 13UI I31j . These moving 
phases have been studied both experimentally jH2] and numerically |33[ I34| 
I35j . Motivated by these results, the dynamical regimes of the moving VL in 
the presence of periodic pinning has also become a subject of interest |T%irHfl 
E3EDE21l2SllSlESll2niEZl- At zero temperature, the dynamical phases of 
vortices driven by an external current with a periodic array of pinning sites 
has been studied in very detail by F. Nori and coworkers |18| ll9j. A complex 
variety of regimes has been found, particularly for N v > N p where the motion 
of interstitial vortices leads to several interesting dynamical phases. 

Moreover, most of the effects of periodic pinning that have been studied 
are related to conmensurability phenomena and the breaking of translational 
symmetry in these systems. Less studied is the effect of the breaking of ro- 
tational symmetry in periodic pinning potentials, in particular regarding 
transport properties. One question of interest is how the motion of vor- 
tices changes when the direction of the driving current is varied. If there is 
rotational symmetry, the vortex motion and voltage response should be in- 
sensitive to the choice of the direction of the current. However, it is clear that 
in a periodic pinning potential the dynamics may depend on the direction of 
the current. For example, in square JJA it has been found experimentally 
that the existence of fractional giant Shapiro steps (FGGS) depends on the 
orientation of the current bias. When the JJA is driven in the [11] direction 
the FGGS are absent, while they are very large when the drive is in the [10] 
direction [2B]- But may be the another example of more recent interest is the 
phenomenon of transverse critical current in superconductors with pinning, 
already mentioned before [201 E01 EH] • Eventhough many numerical results 
have shown the existence of a critical transverse current ^Bl EB HHj there 
has not been experimental measurements of the existence of transverse volt- 
ages and transverse pinning effects before our initial experimental-numerical 
work pij . 

In particular, a very interesting case of conmesurability effect in JJA, 
where the non-equilibrium vortex dynamics could be study are the fully 
frustrated JJA. In the presence of a magnetic field such that there is a half 
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flux quantum per plaquette, / = Ha 2 /&o = 1/2, the JJA corresponds to 
the fully frustrated XY (FFXY) model [201 E2 EBl EH] • The ground state 
is a "checkerboard" vortex lattice, in which a vortex sits in every other site 
of an square grid [HE!- There are two types of competing order and broken 
symmetries: the discrete Zi symmetry of the ground state of the vortex 
lattice, with an associated chiral (Ising-like) order parameter, and the con- 
tinuous £7(1) symmetry associated with superconducting phase coherence. 
The critical behavior of this system has been the subject of several exper- 
imental [SB] and theoretical [SOI EZ1 EHI EDI SD S21 SHI E3| studies. There 
are a Z<i transition (Ising-like) and a f/(l) transition (Kosterlitz-Thouless- 
like) with critical temperatures Tz 2 > -^1/(1) • There is a controversy about 
these temperatures being extremely close |42| or equal [411 1431 IH| . In the 
light of this, it is worth studying the possibility of non-equilibrium Z2 and 
U(V) transitions at large driving currents. Also, the dynamical transitions in 
driven systems studied up to now involve continous (translational or gauge) 
symmetries, and therefore it is interesting to study a system with a discrete 
symmetry. 

In this review we collect our previous works |23| I24[ 1251 1261 127j on the 
study of the dynamical regimes of a moving VL in the periodic pinning of 
a Josephson junction array (JJA) of L x L junctions, at finite temperatures 
for different cases of submatching field (/ = 1/25, / = 1/L 2 and / = 1/2). 
For the case of diluted vortex lattices (/ = 1/25) we obtain a phase dia- 
gram as a function of the driving current / and temperature T. We find 
that when the VL is driven by a low current, the depinning and melting 
transitions can become separated even for a field for which they coincide 
in equilibrium. Moreover, we can distinguish between the depinning of the 
VL in the direction of the current drive, and the transverse depinning in 
the direction perpendicular to the drive. This later case corresponds to 
the vanishing of the transverse critical current in a moving VL at a given 
temperature Tt r , or equivalently, to the vanishing of the transverse supercon- 
ducting coherence. We obtain three distinct regimes at low temperatures: 
(i) Pinned vortex lattice: for < T < T p (I) there is an ordered VL which has 
crystaline long-range order, superconducting coherence (i.e., a finite helicity 
modulus) and zero resistance both in the longitudinal and transverse direc- 
tions, (ii) Transversely pinned vortex lattice: for T P {I) < T < T tr (I) there 
is a moving VL which has anisotropic Bragg peaks, quasi-long range order, 
transverse superconducting coherence and zero transverse resistivity. There 
is a finite transverse critical current. This regime also has strong orienta- 
tional pinning effects in the [1,0] and [0,1] lattice directions, (iii) Floating 
vortex lattice: for 2t r (J) < T < Tm(I) there is a moving VL which is un- 
pinned in both directions and it has quasi-long range crystalline order with 
a strong anisotropy. After our first work in Ref. |2S] , further studies of ther- 
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mal effects in a moving vortex lattice with a periodic pinning array have 
been reported for N v = N p ^ and for N v > N p (213 E2] • Some of these 
results are similar to ours. 

Regarding the breaking of rotational invariance in square J J A, case in 
which the discrete lattice of Josephson junctions induces a periodic egg- 
carton potential for the motion of vortices flj, we will show how the voltage 
response depends on the angle of the current with respect to the lattice di- 
rections of the square JJA. We will show numerically that there are preferred 
directions for vortex motion for which there is orientational pinning. These 
results are in good agreement with experimental results. This leads to an 
anomalous transverse voltage when vortices are driven in directions different 
from the symmetry directions. An analogous effect of a transverse voltage 
due to the guided motion of vortices has been observed experimentally in 
YBCO superconductors with twin boundaries |45j. 

As we mentioned before, in a first stage of our research we have found 
dynamical transitions of the vortex lattice in a JJA with a field density 
of / = 1/25 |2ol 125] : for large currents / there is a melting transition 
of the moving vortex lattice at a temperature higher than the transverse 
superconducting transition: Tm(I) > 2~[/(i)(-0- Interestingly, in the fully 
frustrated case (/ = 1/2) we find that the opposite case occurs in the driven 
FFXY: the order of the "checkerboard" vortex lattice is destroyed at a much 
lower temperature than the transverse superconducting coherence, Tz 2 (I) < 

The remainder of our review is organized as follows. In Sec. 2 we in- 
troduce the theoretical model used for the dynamics of the JJA. We also 
discuss how this model can be mapped to a superconducting film with a pe- 
riodic array of holes. In Sec. 3 we discuss the results for the case of diluted 
vortex lattices dynamics in JJA including the orientational pinning results. 
In Sec. 4 we present the results on the non-equilibrium dynamical regimes 
in fully frustrated JJA. In Sec. 5 we present a discussion comparing the two 
previous sections and the general conclusions. We also provide a detailed 
definition of the adequate periodic boundary conditions for a JJA with an 
external magnetic field and an external driving current in Appendix A, as 
well as the algorithm used for the numerical simulation in Appendix B. 

1.2 MODEL 

1.2.1 Resistively Shunted Junction dynamics 

We study the dynamics of square Josephson junction array (JJA) with LxL 
superconducting nodes (see Fig |l.l|) using the Resistively Shunted Junction 
(RSJ) model for the junctions |46[ I47[ 14*8 } I49[ I50| . The nodes are in the lattice 
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Figure 1.1: Experimental samples: square proximity-effect Pb/Cu/Pb 
Josephson Junction Arrays SEM's image, with an array lattice space a = 
lQfxm and distance between superconducting islands d = l[im (the authors 
thank to Dr. H. Pastoriza and collaborators, Centro Atomico Bariloche, 
Argentina, for this picture.) 



sites n = (n x ,n y ) and their superconducting phases are 9(n). Basically we 
have an XY model, therefore its Hamiltonian is the following: 

U = ~ E ^ cos ^ n + ^ - 0(n) - A,(n)\ . (1.1) 

Using this model, the current flowing in the junction between two su- 
perconducting nodes in the JJA is modeled as the sum of the Josephson 
supercurrent and the normal current: 

J„(n) = Josin^(n) + + ^(„,t) (1.2) 

being Iq the critical current of the junction between the sites n and n + fi, 
(fi = x, y), Rn the normal state resistance and 

0„(n) = 0(n + //) - 0(n) - A^n) = A^(n) - A^n) (1.3) 

the gauge invariant phase difference with 

2vr f( n +f J -) a 



A fl (n) = — A-dl. (1.4) 

*0 Jna 

The thermal noise fluctuations rj^ have correlations 

2kT 

fa / ,(M) V (n',t , )> = —5^,5 n ^5(t - t'). (1.5) 

In the presence of an external magnetic field H we have 

A M xA M (n) = A x (n) - A^n + y) + A,(n + x) - A y (n) 

= 2tt/, (1.6) 



8 



V.I. Marconi, D. Dominguez 



f = Ha 2 /&o and a is the array lattice spacing. If the magnetic field is such 
that / = 1/L 2 it corresponds to have one vortex in the sample, if / = 1/25 
we have an example of diluted vortex lattices, if / = 1/2, there is half flux 
quantum per plaquette and the JJA correspond to the fully frustrated XY 
model. 

We take periodic boundary conditions (PBC) in both directions in the 
presence of an external current I = Iy |35| (See Appendix A). The vector 
potential is taken as 

A,(n,t) = A^n) - a^t) (1.7) 

where in the Landau gauge A x (n) = —2irfn y , Ay(n) = and a^(i) allows 
for total voltage fluctuations under periodic boundary conditions. In this 
gauge the PBC for the phases are [231 HH1 ESI : 

9(n x + L,n y ) = 6(n x ,n y ) 

8(n x ,n y + L) = 9(n x , n y ) - 2ixfLn x . (1.8) 
We also consider local conservation of current, 



(1.9) 



After Eqs. (jl.2ll.9j) we obtain the following equations for the phases (2l?| l23 1 



E3- 



where 



h A2 d9(n) 



2eR N * dt 



-A M .[5 M (n) + 7/^(11,*)] 



(1.10) 
(1.11) 



S^n) = / o sin[A M 0(n) - Aj(n) - a„] , 
and the discrete Laplacian is 

Aj/(n) = /(n + x) + /(n-x) + /(n + y) + /(n-y) 

-4/(n). (1.12) 
The Laplacian can be inverted with the square lattice Green's function G n n / : 



A);G nin , =<y n , n /. (1.13) 
Since we take PBC (see Appendix A), the total current has to be fixed 



by: 



Ix ~ L2 



1 

1? 



sin 6 X (n) + r\ x (n, t) 



sin 6^(11) +r) y (n,t) 



+ 



+ 



h da x 
2eR N ~dT 



(1.14) 



h da 



2eR N dt ' 
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These equations determine the dynamics of a^{t) [33] • For the case of 
a current flowing in the y direction we take I x = and I y = I and for the 
case in which we apply a current / at an angle (j> with respect to the [10] 
lattice direction, I x = I cos cj> and I y = I sin (j>. After Eqs. 1)1. 1011. 1311.14(1 
we obtain the following set of dynamical equations [231 ESI E2] , 



^ = -^G^A^.^nO + r^n',*)], (1.15) 

n' 

= ^-^ w+ ^ (n ' t) ' (L16) 

n 

where we have normalized currents by Jo, time by tj = 2eRNlo/h, and 
temperature by Jo$o/27rfcs- 



1.2.2 Comparison with thin film with a periodic array of 
holes 

Let us consider a superconducting thin film with a square array of holes, 
which act as pinning sites for vortices. There are N p = I? pinning sites 
separated by a distance a. The current density in the superconducting film 
is given by the sum of the supercurrent and the normal current: 



J = 3s + Jn 

j = ^d*-(d*)**] + ^(w-^aV 

m* 1 v ' J 2vrc dt V $o / 

(1.17) 

with D = V + ^^A, VP(r) = |^(r)| exp[i#(r)] the superconducting order 
parameter and a the normal state conductivity. These equations are valid 
everywhere in the film except in the hole regions. If the number of vortices 
N v = BL 2 a 2 is much smaller than the number of pinning sites N p , all 
vortices will be centered in the holes in equilibrium. In this case we can as- 
sume that |^(n)| ~ |^o I is homogeneous in the superconducting film. There- 
fore the dynamics is given by the superconducting phase 0(r), corresponding 
to a London model in a sample with holes. After considering current conser- 
vation V • J = 0, we obtain the London dynamical equations for the phases 
in this multiply connected geometry. Since N v <C N p , we make the approx- 
imation of solving the equations in a discrete grid of spacing a. This means 
that we take as the relevant dynamical variables the phases 0(r n ) defined 
in the sites which are dual to the pinning sites. They represent the aver- 
age superconducting phase in each superconducting square defined by four 
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pinning sites. Therefore, we take the discretization r n = (n x a, n y a) = an. 
[Pinning sites at centered at positions r p = {n x + 1/2,% + l/2)a]. The 
derivatives in the supercurrent are discretized in a gauge-invariant way as 



After doing this, we obtain an equation analogous to (|1.2[) . Now Ia{ n ) 
has to be interpreted as current density normalized by Jo = 2eh\ i f>o\ 2 /ma = 
$o/(8tt 2 A 2 a), time normalized by r = c/ (47T(tA 2 ), and the fraction of vortices 
is / = N v /N p = Ba 2 /<&o. This leads to a set of dynamical equations of the 
same form as Eqs. (|1.15ll.l(i|) . Therefore, we expect that for / <C 1 the model 
for a J JA also gives a good representation of the physics of a superconducting 
film with a square array of holes (meaning that effects of interstitial vortices 
are neglected for N v <C N p ). In other words, we expect that for a low density 
of vortices the specific shape of the periodic pinning potential (being either 
an egg-carton or an array of holes) will not be physically relevant. 

1.2.3 Quantities calculated and simulation parameters 

The Langevin dynamical equations (|1,15I1.16|) are solved with a second or- 
der Runge-Kutta-Helfand-Greenside algorithm with time step At = (0.05 — 
0.1)tj. The discrete periodic Laplacian is inverted with a fast Fourier + 
tridiagonalization algorithm (see the Appendix B for a detail of the algo- 
rithm used). 

The main physical quantities calculated are the following: 
(i) Transverse superconducting coherence: We obtain the helicity modu- 
lus T x in the direction transverse to the current as 



Whenever we calculate the helicity modulus along x, we enforce strict peri- 
odicity in 9 by fixing a x (t) = (see Appendix A). 

(ii) Transport: We calculate the transport response of the J JA from the 
time average of the total voltage as 



D^{r) -> - tf(n + /i) -e 




(1.18) 




V x = (v x (t)) 

Vy = (Vy(t)) 



(da x (t)/dt) 
(da y (t)/dt) 



(1.19) 



with voltages normalized by RnIq- 
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Figure 1.2: Left: schematic representation of a square Josephson junction 
array, the crosses indicate Josephson junctions (JJ), with a rotated applied 
external current. Right: definition of angles with respect to the x and y axis 
of the array. The vector of applied current I = (I x ,I y ) forms an angle <j) 
respect to the x-axis. The vector of electric field E = (V x , V y ) forms an angle 
9 V with respect to the x-axis and an angle 6t with respect to the current I. 



(iii) Vortex structure: We obtain the vorticity at the plaquette n = 
(n x + 1/2, n y + 1/2) (associated to the site n) as |5T] : 

6(n) = —A/j, x nint[^(n)/27r] (1.20) 

with nint[x] the nearest integer of x. We calculate the average vortex struc- 
ture factor as 




(1.21) 



Analyzing above properties we study JJA under different magnetic fields 
applied, (a) Diluted vortex lattices corresponding to / = 1/25. (b) A single 
vortex in the array, / = 1/L 2 . (c) And fully frustrated JJA, correspond- 
ing to have half flux quantum per plaquette, i.e. f = 1/2. We consider 
square networks of L x L junctions, with L = 8, 16, 24, 32, 48, 64, 128 or 
L = 50, 100, 150. We apply dc currents in different directions: (i) I flowing 
in the y direction or (ii) I at an angle <fi with respect to the [10] lattice 
direction, 

I x = I cos 4> 

Iy = Jsinc/>. (1.22) 

We define the longitudinal voltage as the voltage in the direction of the 
applied current, 

Vi = V x cos <j) + V y sine/), (1.23) 
and the transverse voltage 



Vt = —V x sin <f> + V y cos <f>. 



(1.24) 
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From the voltage response, we define the transverse angle as 



tan6 t = V t /V l 



(1.25) 



and the voltage angle as 



tan 9 V = V y /V a 



XI 



(1.26) 



i.e. 6t = V — 4>, see Fig. 11.21 When the vortices move in the direction 
perpendicular to the current, there is no transverse voltage, therefore 9t = 
and 9 V = eft. The case with a dc current applied in the y direction, I = Iy 
{(p = vr/2), will be the more extensively used in this review. 



1.3 DILUTED VORTEX LATTICE DYNAMICS 



In this section we show the study of JJA with a low magnetic field corre- 
sponding to / = 1/25, diluted vortex lattices, for different system sizes of 
L x L junctions, with L = 50, 100, 150. Most of the results are for L = 100, 
except when it is explicitly specified, and for Nt = 10 5 iterations after a 
transient of Nt/2 iterations. 

1.3.1 Transition near equilibrium 

The ground state vortex configuration for / = 1/25 is a tilted square-like 
vortex lattice (VL) j^H, see Fig |1.3f a). We find that this state is stable 
for low currents and low temperatures (in fact, the structure of Fig ll.3f a) 
corresponds to I = 0.01 and T = 0.01). The lattice is oriented in the [4a, 3a] 
direction and commensurated with the underlying periodic pinning potential 
of the square JJA. The structure factor S(k) has the corresponding Bragg 
peaks at wavevectors G in the reciprocal space, as can be seen in Fig ll.3f b). 
When the temperature is increased, the VL tends to disorder and above the 
melting temperature Tm a random vortex array with a liquid-like structure 
factor is obtained, Figs ll.3r c) and inset. 

We find a single equilibrium phase transition (/ = 0) at Tm ~ 0.050 db 
0.003, which is in agreement with the melting temperature obtained by Franz 
and TeitelHH and Hattel and WheatleyUHl for / > 1/30. 

We now apply a very low current, I = 0.01, in order to study the 
near-equilibrium transport response simultaneously with other quantities 
like structure factor and helicity modulus. We find a phase transition at a 
temperature Tm(I) ~ 0.046 ± 0.001, which is slightly lower than the equi- 
librium transition. In Fig ll.4f a) we see that there is a large jump in the 
resistance R = V/I at T c , in good agreement with the first-order nature of 
the equilibrium transition ^3] . The onset of resistivity is a signature of a 
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Figure 1.3: Vortex configuration for a low current I = 0.01 <C ^ c (0) and a 
low temperature, T = 0.01 < T/^: (a) tilted square vortex lattice, oriented in 
the [4a, 3a] direction and (b) corresponding structure factor S(k). At high 
temperature, T = 0.05 > Tm- (c) disordered vortex array. Inset: vortex 
liquiddike structure factor. 
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Figure 1.5: Zero temperature curve. There is a critical current I c (0) = 
0.114 ±0.002. 
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Figure 1.6: TV curves for different temperatures. Linear scale, crossing 
point at I* = 0.165. Inset: Log-log plot of R = V/I vs. I for the same 
temperatures. 



depinning transition in the direction of the drive. This occurs simulta- 
neously with a melting of the vortex lattice, corresponding to the vanishing 
of Bragg peaks, as shown in Fig ll.4f b) for the two first reciprocal lattice 
vectors d = ^(-4/25,-3/25) and G 2 = ^(-3/25,4/25). In the direc- 
tion of the current drive the helicity modulus is ill-defined since total phase 
fluctuations are allowed (see Appendix A). However, in the perpendicular 
direction to the drive the helicity modulus can be calculated, and it is 
a measure of the transverse superconducting coherence. As we can see in 
Fig ll.4f c). transverse superconductivity also vanishes at Tm{I)- Above Tm, 
we find that the T X (T) has large fluctuations around zero. 



1.3.2 Transport properties 

Let us now study the transport properties for larger currents. We calculate 
the current-voltage (IV) characteristics for different temperatures as well as 
the dc resistance R = V/I as a function of temperature for finite currents. 
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Figure 1.7: < V y > /I y vs temperature curves for different dc currents (Iy). 



The zero temperature IV curve has a critical current, I c (0) = 0.114 ± 
0.002, see Fig ll.5| which corresponds to the single vortex depinning current 
in square JJA [I], with the typical square root depence at the onset. Similar 
behavior has been reported for zero temperature IV curves for low values 
of / [531 154] . Above I c (0) there is an almost linear increase of voltage, 
corresponding to a "flux flow" regime, where there is a fastly moving VL. The 
structure factor of the T = moving VL is the same as the corresponding one 
of the pinned VL [Fig ll.Sf a)]. The presence of periodic boundary conditions 
in our case prevents the occurrence of random or chaotic vortex motion 
near the critical current, as reported in early simulations with free boundary 
conditions |461 147] . In what follows we will restrict our analysis for currents 
I < 0.4, where the collective behavior of the VL is the dominant physics 
(at I ~ 1 there is a sharp increase of voltage when all the junctions become 
normal, and V ~ RnI for I 1). 

The IV curves for finite temperatures are shown in Fig ll.61 For temper- 
atures below Tjvf there is a nonlinear sharp rise in voltage which defines the 
apparent critical current I C (T). For example, we can obtain this I C (T) with 
a voltage criterion, which we choose as V < N j At = 10~ 4 . In this case, we 
find that I C (T) decreases with T, vanishing at Tm- It is interesting to point 
out that all the IV curves for different temperatures have a crossing point 
at I* = 0.165, see Fig ll.6l A crossing in the IVs has also been reported in 
experiments in amorphous thin films [55] . For temperatures T > Tm the IV 
curves tend to linear resistivity for low currents. This is shown in the inset 
of Fig ll.6l in a log-log plot of R(I) = V/I vs /, where we see that R(I) tends 
to a low current finite value for T > Tm, while it has a strong nonlinear 
decrease for T < Tm- 
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How to get "schematically" the transverse temperature depin- 



Let us now study the dc resistance R = V/I as a function of temperature 
for a given applied current in the y-direction (Fig ll.7jl . We start with the 
perfectly ordered VL as an initial condition at T = and then we slowly 
increase the temperature, keeping / constant. For currents below the T = 
critical current, / < I c (0), the dc resistance is negligibly small at low T, 
and it has a steep increase at a depinning temperature T p (I), corresponding 
to the onset of vortex motion. The depinning temperature decreases for 
increasing currents, and the values of T P (I) are coincident with the apparent 
critical currents I C (T) obtained from the IV curves. For currents higher 
than ic(0), there is always a large and finite voltage for any temperature. 
If / c (0) < I < I* , the R(I) increases slightly with T tending to a constant 
value for large T, while for I > I* the R(I) decreases with T. 

1.3.3 Transverse depinning 

What is its response to a small current in the transverse direction when 
the driven vortex lattice is moving? Is the vortex lattice still pinned in the 
transverse direction? Is there a transverse critical current for a moving VL? 
The idea of a transverse depinning current was introduced by Giamarchi 
and Le Doussal in Ref. 3Q] for moving vortex systems with random pinning 
at zero temperature. The possibility of such a critical current was later 
questioned by Balents, Marchetti and Radzihovsky where it was shown 
that this is not true for any finite temperature in random pinning] however 
a strong nonlinear increase of the transverse voltage was predicted at an 
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Figure 1.9: Transverse depinnig temperature determinations. With a small 
transverse current applied, I x = 0.01, and I y = 0.16, < Vtr/hr > y s T. 
Inset: Tr — IV curve and < V y > vs Itr/Iy at T = and high longitudinal 
current, I y = 0.16. 



"effective" transverse critical current. In the case of periodic pinning it is 
more clear that a transverse critical current will exist at T = since it is 
a commensurability effect. This has been found in the T = simulation 
work of Reichhardt et al ^Sl • It is also possible that this transverse critical 
current will still be non-zero at T ^ in periodic pinning. In fact, we have 
found in our previous work j2H] that there is a thermal transverse depinning 
in a periodic system, and we will now analyze this behavior in detail. 

First, a high longitudinal current / > I c (0) is applied at zero tempera- 
ture. Then, a current It r is applied in the transverse direction (see Fig |1.8l 
scheme). In this way, a transverse current-voltage characteristics can be 
obtained for each /. This is shown in the inset of Fig |1.9l for / = 0.16. We 
clearly see that there is a finite transverse critical current I Ct t r ~ 0.12, which 
is of the order of the single vortex pinning barrier. We also show in the inset 
of Fig ll.91 how the longitudinal voltage V changes when the transverse cur- 
rent is varied. For It r < I c ,tr the longitudinal voltage V is almost constant. 
At I Ct tr there is a fast decay of V. When I tr = I we also have V = Vt r 
as expected for a drive at a degree of tt/4. Later, for I tr > I the vortex 
lattice becomes pinned in the other direction, since now the directions of 
"longitudinal" and "transverse" current are interchanged. 

Another possible measurement is to study thermal transverse depinning. 
In this case, we start with a longitudinal current I, then a small transverse 
current is applied, It r <C I c ,tr, and the temperature is slowly increased. In 
this way, we can measure a transverse resistance Rt r = Vtr/hr like it was 
showed schematically in Fig |1.8l In Fig |1.9l we plot this result for / = 0.16 
and Itr = 0.01. We find that for finite low temperatures Rt r is negligibly 
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Figure 1.10: Transverse IV curves for different temperatures: (a) low cur- 
rent Iy < I c (0), I y = 0.06, (b) high current I y > J c (0), I y = 0.16. 



small within our numerical accuracy and it has a clear onset at a transverse 
critical temperature Tt r . 

Let us now see how these results depend on the longitudinal current / 
and temperature T. We have calculated the transverse IV (Tr-IV) curves 
for different I and T. In Fig lLlOf a) we show the Tr-IV curves for I = 0.06 
(low current regime) and in Fig ll.lOf b) for / = 0.16 (high current regime). 
In both cases, there is a clear change of behavior in the Tr-IV curves when 
going through a characteristic transverse critical temperature Tt r (I). For 
low temperatures there is a transverse critical current which tends to vanish 
when T aproachs Tt r (I) from below. In contrast, for T > Tt r there is a linear 
resistivity behaviour. 

The transverse resistivity Rt r = Vt r /Itr as a function of temperature 
was calculated for different longitudinal currents fFig ll.l In all cases 
there is an onset of transverse response at a given temperature Tt r {I). At 
low currents, I < I c (0), the transverse depinning temperatures are almost 
constant, Tt r ~ 0.02 tending to increase slowly with /, see Fig ll.lll On the 
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Figure 1.11: Vtr/hr vs T curves for different dc currents applied in the y 
direction and small transverse current applied, I x = 0.01: for low current 
I y < J c (0). Inset: for high current I y > I c (0). 



other hand, for / > I c (0) the transverse depinning temperatures increase 
clearly with /, see the inset in Fig Jl.llI 

1.3.4 Non-equilibrium regimes 

We will now study the different non-equilibrium regimes of vortex driven 
lattices and characterize their possible dynamical transitions. The approach 
we will follow in this subsection is to have a fixed current applied in the 
system and vary the temperature. In this way, we look for the possible 
transitions as a function of T in a similar way as was done near equilibrium 
in Sec. 1.3.1. 

A few similar studies were done previously in related systems. In Ref. [HI] 
the melting transition of a moving vortex lattice in the three dimensional 
XY model was studied in this way. In this work a first order transition 
was found as a function of temperature in a strongly driven vortex lattice. 
In Ref. j5U] a current driven two dimensional JJA at zero field was studied. 
The possibility of a transition as a function of temperature for finite currents 
below the critical current was analyzed in this case. 

In the following, we will separate our study in three ranges of current: 
low currents, I < 0.04 , intermediate currents 0.04 < I < I c (0) , and high 
currents I > / c (0) . 

Low currents 

We show the results for a low current in Fig ll.l2l for / = 0.03. The longitu- 
dinal dc response, V/I, is negligibly small at low T and later it has a sharp 
increase of two orders of magnitude: this defines the depinning temperature, 
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Figure 1.12: For low current I < J c (0), / = 0.03: (a) < V y > /I y vs. T, (b) 
5(Gi) (o) and 5(G2) (*) vs. T, inset: size effect in S(Gi), (c) T x vs. T. 



T p [Fig ll,12T a)]. At higher temperatures, T > T p , the resistance is weakly 
T-dependent. 

Below T p the vortex lattice is pinned and its structure is similar to the 
ground state: a vortex lattice commensurate with the underlying square 
array and tilted in the [4a, 3a] direction. Above T p the vortex lattice is 
moving and it has an anisotropic structural order. If we analyze the structure 
factor in two different reciprocal lattice directions S(Gi) and S(G2), we can 
see this clearly (Fig ll.l2f b)). Below T p the VL structure is isotropic and 
S(Gi) = S(G2)- Right above T p the height of the peaks decreases with 
temperature, and the structure of the depinned VL is clearly anisotropic, 
S(Gi) 7^ S(G2)- Finally, at a melting temperature Tm the peaks vanish, and 
the vortex lattice melts into a vortex liquid. In Fig ll.^f b) inset we compare 
the behavior of the Bragg peaks for two system sizes L = 50, 100. We see 
that the S(G\ t 2) are size independent below T p as it should be expected for 
a pinned phase [T4*] . 

The helicity modulus in the direction perpendicular to the current, T x , 
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Figure 1.13: For I < I c (0), / = 0.06: (a) < V tr > /I tr vs. T, (b) < V y > / I y 
vs. T, (c) 5(Gi) (o) and 5(G 2 ) (★) vs. T, (d) T x vs. T. 



decreases very slowly for T < T p . Above T p , T x has a faster decay with 
important fluctuations and tends to vanish at Tm ■ For T > Tm, T x oscillates 
around zero. Therefore, for a small finite current the depinning and melting 
transitions become separated with T p < Tm- 

Intermediate currents 

At intermediate currents, 0.04 < I < I c (0), a new transition appears: the 
transverse depinning of the moving vortex lattice. As discussed in Sec. 1.3.3, 
one can measure transverse depinning by applying a small transverse current 
while the VL is driven with a fixed longitudinal current. This is shown in Fig. 
11.13( a) for the case of I = 0.06 and a small transverse current, It r = 0.01. 
We see that there is an onset of transverse voltage at Tt r = 0.019. We can 
also see that this transverse depinning temperature is above the depinning 
temperature T p for longitudinal resistance [Fig ll. 13( b). T p = 0.013], and 
below the melting temperature Tm for the vanishing of the Bragg peaks 
[Fig ll. 13( c). Tm = 0.037]. Therefore, this transition occurs at an inter- 
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Figure 1.14: For intermediate current I < I c (0), I = 0.06: size effect in 
5(d). 



mediate temperature between the depinning and the melting transitions, 
T p < T tr < TV. We also show in Fig ll.l.Sf d) that the helicity modulus 
begins to fall down slowly at Tp r , while for T% r < T < Tm it has strong 
fluctuations, being difficult to interpret its behavior in this case. 

One can see that the intensity of the Bragg peaks has a greater depen- 
dence with size for T tr < T < Tm when compared with the T p < T < T tr 
regime in Fig ll.141 For large temperatures T > Tm, in the liquid phase, the 
value of S(Gi t 2) in general is strongly size dependent, since it should go as 
L~ 2 |T1]. On the other hand, for T p < T < Tm we find that the intensity 
of the Bragg peaks is weakly dependent on system size. The clear change 
of behavior of the size dependence gives a good criterion to determine Tm 
(FigHHI). 

It is interesting to study in more detail the behavior of the structure 
factor through all these transitions. In Fig |1.15l we show examples of the 
structure factor S(k) at temperatures in the different regimes. In the pinned 
phase, the S(k) is nearly the same as in the ground state with delta-like 
Bragg peaks, see Fig ll,15f a). For the moving VL, we can see that there is 
less anisotropy in the transversely pinned regime T p < T < Tt r [Fig |l.l5T b)] 
than in the floating regime Tf r < T < Tm [Fig |1.15T c-d)]. Moreover, near 
Tm the VL structure becomes strongly anisotropic, with the peak at Gi 
much larger than the peak at G2 (see Fig lLlSt d) and Fig |1.14|) . 

The anisotropy of the Bragg peaks of the moving VL (in the regimes at 
T p < T < Tm) has two characteristics: (i) the width of the peaks increases 
with T in the direction of the applied current (the direction perpendicular to 
the vortex motion), and (ii) the height of the peaks decreases in the direction 
of vortex motion. This can be observed in the sequence of structure factors 
shown in Fig |1.15l 
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Figure 1.15: Intensity plot of the structure factor 5(G) at I < I c (0), I = 0.06 
for different temperatures: (a) T < T p , T = 0.006, (b) T p < T < T tr , 
T = 0.0185, (c) <T < T M , T = 0.024, (d) T < T M , T = 0.035. 
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Figure 1.16: For high current I > I c (0), / = 0.16: (a) < V tr > /I tr vs. T, 
(b) 5(Gi) (o) and 5(G 2 ) (*) vs. T, (c) T x vs. T. 



High currents 

In the case of currents larger than I c (0), the vortex lattice is already de- 
pinned at T = 0. As we have seen in Sec. 1.3.3, this zero-temperature 
moving VL has a finite transverse critical current, and therefore it is pinned 
in the transverse direction. When we slowly increase temperature from this 
state, we find that the transverse resistive response is negligible for finite low 
temperatures. At a temperature Tf r there is a jump to a finite transverse 
resistance Rt r = Vtr /Itr- For example, this is shown for I = 0.16 with a 
small transverse current, I x = 0.01 in Fig |1.16f a). The vortex lattice has an 
anisotropic structural order for all temperatures, ie, S[G\) ^ 5(G2), and 
the height of the Bragg peaks vanishes at Tm, Fig |1.16f b). The transverse 
helicity modulus (T x ) is almost constant for low temperatures and starts to 
decrease at T tr presenting strong fluctuations for T > T tr , Fi jl.lHf c). In 
Fig |1.17l we show 5(G2) for different sizes, L = 50, 100, and we see that Tm 
is size independent. Similar behaviour is found for S'(Gi). 

The analysis of the heights of the peaks in 5"(k) as a function of system 
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Figure 1.17: For high current / > 7 C (0), / = 0.16: size effect in S(Gz) vs. 
T. 



size L is a good indicator of the translational correlations in the system. This 
dependence is well known for two dimensional lattices ^1] . For a pinned solid 
S(G) ~ 1, for a floating solid, S(G) ~ L~^{I,T) with < m (I,T) < 2, 
being this dependence a signature of quasi-long range order, and for a normal 
liquid S(G) ~ L~ 2 . To assure the existence of algebraic translational 
correlations, we have done this scaling study for currents 0.02 < I < 0.2 and 
different temperatures for system sizes of L = 50, 100, 150. For the cases 
corresponding to the pinned regime (T < T p ) we found t\g ~ 0, as expected. 
In the transversely pinned regime, we show a case in Fig ll.lSf a). we find a 
power law fitting with very small values of t]g{I, T). In the floating regime, 
we find larger values of t]g, we show a case in Fig ll.lSf b). In all the cases we 
have obtained that < t/g{I, T) < 2 for T p < T < Tm- Therefore, this finite 
size analysis shows the existence of quasi-long range order in the moving VL. 
Also, we find that ijg 1 > rjG 2 f° r an currents and temperatures. The power- 
law exponent t)g can be studied at a constant current, as a function of 
temperature. This is shown in the inset of Fig lLlSf b). for / = 0.16 and 
different reciprocal lattice vectors. The exponent r\ is finite for the complete 
temperature range. For T < Tj r it has a small value t)g ~ 0.01. It has a 
fast increase near Tf r , where there is also a clear difference between tjg 1 and 
r]G 2 ■ Finally, it reaches a value of r/ = 2 near T = Tm ■ 

1.3.5 Orientational pinning effects 

A very interesting characterization of the different regimes showed in the 
precedent subsection can be obtained by studying the effects of varying 
current direction |24l 1561 I57j . We apply a current / at an angle (j) with 
respect to the [10] lattice direction, I x = Icoscj), I y = Isincj) (see Sec. 1.2.3). 
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Figure 1.18: Finite size analysis and power law fit of 5(G) ~ L VG (I, T) for 
a moving vortex lattice at I > I c (0), I = 0.16. We obtain: (a) for T < T tr , 
T = 0.02, r) Gl = 0.023 (*) and n G , 2 = 0.013, (A), (b) for T tr < T < T M , 
T = 0.035, rjGi = 0-471 (★) and r)o 2 = 0.34 (A). Inset: at same I, rj(I, G) 
vs T, for Giand G 2 . 



We study the voltage response when varying the orientation <fi of the drive 
while keeping fixed the amplitude I of the current. In the parametric curves 
of V y ((f>) vs. V x ((f)) can we analyze the breaking of rotational symmetry in 
the different regimes of / and T. In the case of rotational symmetry this kind 
of plot should give a perfect circle. However, the square symmetry of the 
Josephson lattice will show up in the shape of the curves. In what follows 
we start studying the motion of single vortices in JJA (1), to show later 
the breaking of rotational invariance in diluted vortex lattices (2) and how 
this effect could be used to characterized the dynamical regimes described 
previously. The first case of a single vortex is studied before for its simplicity. 
In order to understand the phenomenon of breaking of rotational symmetry 
is better to begin without including the collective effects of vortex lattices. 
In addition, experimental evidences were found simultaneously for Dr. H. 
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Figure 1.19: Simulation results of the parametric curves V y {4>) vs - Vx (<!>)• 
(a) I = 0.2, for (going outwards from the center) T = 0, T = 0.01 and 
T = 0.05; (b) I = 0.6 for T = 0.3, T = 0.5, T = 0.6, T = 1.0 and T = 1.4; 
(c) / = 1.2 for T = 0, T = 0.05, T = 0.1, T = 0.3 and T = 1.0. (d) T = 0.05 
for (going outwards from the center) / = 1.0, I = 1.2, / = 1.4 and / = 1.6; 
(e) T = 0.2 for / = 1.0, / = 1.2, I = 1.4 and / = 1.6; (f) T = 1.0 for I = 0.2, 
I = 0.4, I = 0.6 and I = 0.8. 



Pastoriza and collaborators [21] for this case. And very nice agreement 
between numerical and experimental results was found. 

Breaking of rotational invariance in JJA with a single vortex 

The square lattice has two directions of maximum symmetry: the [10] and 
the [11] directions (and the ones obtained from them by tt/2 rotations), 
which correspond to the directions of reflection symmetry. When the current 
bias is in the [10] direction, the angle of the current is (j) = 0, and we call it 
a "parallel" bias. When the current bias is in the [11] direction, the angle of 
the current is 4> = 7r /4 = 45°, and we call it a "diagonal" bias (see squeme 
in Fig. Oil . 
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In the case of the parallel bias we find that the transverse voltage is zero 
(in agreement with the reflection symmetry), see details in Ref. [21]. This 
corresponds to vortex motion in the direction perpendicular to the current 
(9t = 0). A perfect agreement with experimental results was found [24]. 

In the IV curve for the longitudinal voltage we find numerically a critical 
current corresponding to the single vortex depinning /j. 10 ' = 0.1 [U I24j . 
In the case of the perfect diagonal bias, which is only attainable in the 
numerical model, we obtain similar results as in the parallel bias case. The 
transverse voltage is zero and therefore the vortex moves perpendicular to 
the current (9t = 0). In agreement with the prediction of Halsey |58| . the 
numerical IV curve for V\ has a critical current of ij. 11 ^ = ^I [ c 0] = 0.1414. 
The onset of the resistive regime is also multiplied by a factor of \/2. 

For orientations different than the symmetry directions, we always find a 
finite transverse voltage (as well as was found in experiments |2U)- ^ n order 
to see this, we study the voltage response when varying the orientation 
4> of the drive while keeping fixed the amplitude I of the current. If we 
analyze the transverse angle 6t = arctan(Vj/V^) as a function of the angle of 
the current (f> we observe that 0t vanishes only in the maximum symmetry 
directions corresponding to angles (j) = 0, ±45°, ±90°, .... Furthermore, we 
see that for orientations near <f> = 0, the transverse angle basically follows the 
current angle: 9t ~ — 4>- This is an indication that vortex motion is pinned 
in the lattice direction [10], since V y m 0, meaning that the voltage angle 
is 8 V ~ 0. Whenever the voltage response is insensitive to small changes in 
the orientation of the current, we will call this phenomenon orientational 
pinning. On the other hand, near <p = 45° the transverse angle changes 
rapidly. 

A more direct evidence of the breaking of rotational symmetry can be 
seen in the parametric curves of V y {4>) vs - V x (4>). In Fig. 11.191 we plot the 
values obtained for the voltages V y and V x when varying the orientational 
angle 4> for different values of the current amplitude I and the temperature. 
In the case of rotational symmetry we should have a perfect circle. In the set 
of plots of Fig. ll.l9r a-c). the current amplitude is fixed and the temperature 
is varied. In Fig. Il,19f a) we have I = 0.2, near the onset of single vortex 
motion in the regime. In this case most of the points are either on the 
axis V x = or on the axis V y = 0, indicating strong orientational pinning 
in the lattice directions [10] or [01]. When increasing T the orientational 
pinning decreases and the length of the "horns" in the x and y axis decreases. 
Fig. 11.19T b) corresponds to I = 0.6 , when the vortex is moving fast. In this 
case the horns have disappeared and orientational pinning is lost. However, 
the breaking of rotational symmetry is still present in the star-shaped curves 
that we find at low T. The dip at 45° in the stars are because in this 
direction the voltages are minimum, since the critical current is maximum 
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in this case, 4 1 = 0.1V2. When increasing the temperature, the stars 
tend to the circular shape of rotational invariance. Above the onset of the 
resistive regime the "horned" curves reappear. In this case the orientational 
pinning corresponds to the locking of ohmic dissipation in the junctions in 
one of the lattice directions, either [10] or [01]. Once again, when increasing 
T the horned structure shrinks, and the curves evolve continuously from 
square shapes to circular shapes. 

The variation with current of the rotational parametric curves for a fixed 
temperature is shown in the numerical results of Figs. I1.19f d-f). At a low 
temperature, T = 0.05 we clearly see the horned structure of the curves for 
almost all the currents and even for large currents the circular curves have 
"horns", see Fig. ll,19T d). At an intermediate temperature, T = 0.2, there 
are still some signatures of the orientational pinning [Fig. ll.lST e)]. while for 
T = 1.2 all the curves are smooth and rounded with a slightly square shape 

[Fig. Egp)]. 

Transverse voltage near the [1,1] direction 

In Fig. I1.20f a) we show experimental voltage-current characteristics for 
an array of 100 x 1000 junctions at a low temperature T = 1.25K and at a low 
magnetic field (results from Dr. Hernan Pastoriza's Group, Argentina). The 
current is applied nominally in the [11] direction, but a small misalignment 
is possible in the setup of electrical contacts, therefore <f> = 45° ± 5°. We 
see that for low currents there is a very large value of the transverse voltage 
Vt, which is nearly of the same magnitude as the longitudinal voltage V[. 
The transverse voltage is maximum at a characteristic current I m . Above 
I m , Vt decreases with increasing current while V\ increases. It is remarkable 
that these results are very different from the IV curve obtained numerically, 
where Vt = at <fi = 45°. However, if we assume a misalignment of a few 
degrees with respect to the [11] direction we can reproduce the experimental 
results. In Fig. ll^Of b'l we show the IV curves obtained numerically for 
(j> = 40° and T = 0.02. We see that for low currents Vt is close to Vf 
Vt < Vi, similar to the experiment, and later Vt has a maximum at a current 
J m ~ 1/ cos (/) ~ 1.3. This corresponds to the current for which the junctions 
in the x-direction become critical (I x = 1). The highest transverse voltage 
can be obtained for orientations near <j> = 45° as we explained before. 

Therefore a slight misalignment of the array from the [11] direction is 
useful for studying both experimentally and numerically the behavior of the 
transverse voltage as a function of current and temperature. 

Transverse voltage near the [1,1] direction in disordered JJA 
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Figure 1.20: (a) Experimental results of longitudinal and transverse volt- 
age for a current near the [11] direction, <p = 45° ± 5 at T = 1.25 K. (b) 
Longitudinal and transverse voltage obtained numerically for cj) = 40° at 
T = 0.02hIo/2ekB- (c) Idem (b) for an extended current range. 



Before we studied the breaking of rotational invariance in square JJA|24| 
at finite temperatures analyzing transport properties and we obtained that 
the transverse voltage vanishes only in the directions of maximum symmetry 
of the square lattice: the [10] and [01] direction (parallel bias) and the [11] 
direction (diagonal bias). For diagonal bias this result is highly unstable 
against small variations of the angle of the applied current, leading to a rapid 
change from zero transverse voltage to a large transverse voltage within a 
few degrees. Now we will show that a small amount of disorder induces 
finite transverse voltage in square JJA with a perfect diagonal bias. The 
transverse voltage as a function of current presents a peak which does not 
depend strongly with disorder for moderate disorder strength. This result is 
experimentally relevant since samples always could present a small amount of 
disorder due to tiny fabrication defects in addition to a possible misalignment 
in the setup of electrical contacts. 
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Figure 1.21: Numerical IV characteristics in a system of 32 x 32 junctions 
with diagonal bias, at T = 0.02, and 5% of disorder compared with transverse 
voltage without disorder. Inset: transverse voltage peak for 5%, 10% and 
15% of disorder. 




We simulate an square JJA with disorder in the critical currents using 
the RSJ model introduced before in Sec. 1.2.1. The modification to the 
current flowing in the junction between two superconducting islands is the 
following 

J„(n) = ^(n)Josin^(n) + + Vfl (n,t), (1.27) 

where disorder is introduced through the factor <?^(n) = 1 + A(£„(n) — 1/2). 
£ is a uniform distributed random variable in [0, 1]. Therefore Jn(l — y) < 
5^(n)/o < ^o(l+v)- Again we calculate the time average of the total voltage 
in both directions, V x and V y , and transverse and logitudinal voltages from 
them, as in Ea. Oil and HTM 

In Fig |1.2ll we show current- voltage characteristics, at T = 0.02 (in units 
of the Josephson energy, ksEj) obtained numerically for an array of 32 x 32 
junctions without (A = 0) and with disorder(A = 0.1, i.e. ±5% variation) 
in the critical currents. The dc current is applied in the diagonal direction 
((f) = 45°), [1,1]. We clearly see the appearance of finite transverse voltage, 
Vt r 7^ 0, with a 5% of disorder in the critical currents. For low currents 
longitudinal voltage is equal to transverse voltage and they start to sepa- 
rate when the vortex-antivortex pairs density increases considerably upon 
increasing current. Vt r presents a maximum around I = Iq and for larger 
currents tends to zero. This behaviour is in contrast to the one obtained 
before without disorder, where Vt r = for all currents. We calculate IV 
characteristics as a function of the intensity of disorder and we observe, as 
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Figure 1.22: Numerical simulations: IV curves for A = 0.1 (Vi(°) and 

Vt r (A)). 



we show in the inset, that the transverse voltage is almost independent of 
A in the range considered, 0.1 < A < 0.6. 

In Fig ll.22l we show in a log-log scale our previous results in order to 
compair them with the experimental results, obtained in an array of 100 x 
1000 proximity-effect Pb/Cu/Pb junctions, at a low temperature T = 1.25K 
and at a low magnetic field (see Fig. 11.20(1 . Qualitatively we get a good 
agreement between experiments and simulations now in the case of diagonal 
bias, in opposite to the results obtained without disorder and diagonal bias. 

The simulation results presented up to this point were focused in the 
motion of a single vortex in the periodic pinning of a square JJA. The 
vortex collective effects, for fields / > 1/L 2 , will be discussed below. 

Breaking of rotational invariance in diluted vortex lattices 

In the previous subsection we showed that the "diagonal" [11] direction 
is unstable against small changes in the angle 4>, while the [10] and [01] 
directions are the preferred directions for single vortex motion [23]. This 
shows as "horns" in parametric V y vs. V x plots, which are finite segments of 
points lying in the x or the y axis. This implies the existence of transverse 
pinning in these directions (thus, it corresponds to orientational pinning). 
Now we return to the case of diluted vortex lattices (/ = 1/25) and we 
perform the same analysis as before for the different regimes of the moving 
VL obtained in Sec. 1.3.4. 

In Fig |1.23l we plot the voltages V y and V x when varying the orienta- 
tional angle <j) for different current amplitudes / and temperatures T. In 
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Figure 1.23: Parametric curves V y ((j)) vs. V x (4>) for different dc current 
applied and temperatures: (a) \I y \ = 0.06 and T = 0.015, (b) \I y \ = 0.06 
and T = 0.03, (c) \I y \ = 0.16 and T = 0.015, (d) \I y \ = 0.16 and T = 0.035. 



Fig ll.23f a) we have / = 0.06 and T = 0.015, corresponding to the regime of 
a transversely pinned lattice. In this case most of the points are lying either 
on the axis V x = or on the axis V y = 0, indicating strong orientational 
pinning in the symmetry lattice directions [10] and [01]. When increasing T 
the orientational pinning decreases and the length of the "horns" in the x 
and y axis decreases. Fig ll.23r b) shows results for / = 0.06 and T = 0.03, 
which correspond to T > Tt r when there is a finite transverse resistance. In 
this case the horns have disappeared and orientational pinning is lost. How- 
ever, the breaking of rotational symmetry is still present in the star-shaped 
curves. Also in the high current regime, Fig |1.23f c'). for a low temperature 
T = 0.015 (T < Tt r ) we see that there is orientational pinning with the 
presence of horns, which again disappears for T > Tt r as it is shown in 
Fig |1.23f d'l at T = 0.035. Finally for T ^> Tm, deep inside in the liquid 
phase, the stars tend to the circular shape of rotational invariance. 

Therefore we have shown that orientational pinning is a useful phe- 
nomenon to characterize the different dynamical regimes in JJA, both nu- 
merically and experimentally. 
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Figure 1.24: I — T Phase diagram for / = 1/25. Tm{I) line is obtained from 
5(G) vs. T curves (•). T P {I) line is obtained from IV curves (*), from 5(G) 
vs. T curves (•) and from (V y ) vs. T curves (A). Tf r (I) curve is obtained 
from (Vtr) vs. T curves (♦). 



1.3.6 Summary and discussion 

With all the information of the previous sections we obtain the current- 
temperature phase diagram shown in Fig |1.24l For finite currents we have 
been able to identify three different regimes, a pinned VL for T < T P (I), 
a transversely pinned VL for T p (I) < T < Tf r (T), and a floating VL for 
Tt r (I) < T < Tm{L). It is, however, difficult to define if the temperatures 
T p , T tr ,T]\,f correspond to either phase transitions or to dynamical crossovers. 
The comparisons we have made of the behavior of voltages (longitudinal and 
transversal), structure factor and helicity modulus show that something is 
happening at these temperatures. Also, the comparison for different system 
sizes of the behavior of 5(G) suggest transitions for T p ,T tr ,TM- The trans- 
verse helicity modulus has strong fluctuations for Tt r < T < Tm- These 
fluctuations are not reduced when increasing the simulation time in a fac- 
tor of 10. This could mean that actually the range of temperatures of 
Tf r < T < Tm corresponds to a long crossover region towards a liquid state. 
Also, one may question if the use of the helicity modulus for these far from 
equilibrium states is correct, since T x has been defined from the response 
of the equilibrium free energy to a twist in the boundary condition. Most 
likely, the large fluctuations in T tr < T < Tm are due to the very unstable 
and history-dependent steady states we find in this regime. For example 
the steady states obtained when increasing temperature from T = at fixed 
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current differ from the steady states obtained when increasing current at 
fixed temperature in this regime. They differ in the degree of anisotropy 
and intensity of the Bragg peaks in the structure factor. We think that this 
reflects the fact that the VL is unpinned in all directions and the orientation 
of the moving state will depend on the history for this regime. 

In any case, we can clearly distinguish different dynamical states which 
define the regimes shown in the phase diagram of Fig |1.24l From the ex- 
perimental point of view these transitions (or crossovers) are possible to 
measure. The depinning temperature T p can be obtained from resistance 
measurements at finite currents. The transverse current-voltage character- 
istics and the transverse resistivity can be measured for different longitudinal 
currents and temperatures, and therefore T tr could be obtained. The helic- 
ity modulus can be measured with the two-coil technique pj. It could be 
very interesting to see the results of such measurements at finite currents. 
The structure factor and melting transition can not be measured directly. 
However, in the presence of an external rf current, the disappearance of 
Shapiro steps can indicate the melting transition since they are sen- 
sitive to the translational order of the vortex lattice |6Uj . Therefore, we 
expect that the results obtained here could motivate new experiments for 
fractional or submatching fields in Josephson junction arrays as well as in 
other superconductors with periodic pinning. 

After our articles (see Ref . 23l . |24j . |25j ) several experimental works ap- 
peared. First of all, D. Shalom and H. Pastoriza designed a new two coil 
kinetic inductance technique in order to measure transport anisotropy in- 
duced by the applied current itself. They found for first time experimental 
evidence for the anisotropic character of the current driven diluted vortex 
lattices states in JJA [B^. Later they improved the previous technique, 
implementing rectangular coils with high aspect radio which are now litho- 
graphically fabricated on top of the sample, separated for isolation layers 

M- 

Regarding numerical works, it is very interesting to mention that re- 
cently the disappearance of Shapiro steps was used to indicate the melting 
transition for another submatching field, / = 1/2 |63j . system which we are 
going to discuss in next section. 

About our results on orientational pinning we would like to do a series 
of remarks and comments. In the egg-carton potential of a square JJA there 
are pinning barriers for vortex motion in all the directions. The direction 
with the lowest pinning barrier is the [10] direction. Therefore the strong 
orientational pinning we find here is in the direction of the lowest pinning 
for motion, i.e. the direction of easy flow for vortices. The presence of a 
strong orientational pinning leads to a large transverse voltage when the 
systems is driven away from the favorable direction, to the existence of a 
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critical angle and to a transverse critical current. On the other hand, the 
[11] direction is the direction of the largest barrier for vortex motion in the 
egg-carton potential. In this case, the behavior is highly unstable against 
small variations in the angle of the drive, leading to a rapid change from 
zero transverse voltage to a large transverse voltage within a few degrees. 
Any misalignment of the current/voltage contacts as well as disorder in the 
junctions critical currents [27j can lead to a large transverse voltage for 
arrays with a diagonal bias. This explains the transverse voltage observed 
experimentally in JJA driven near the diagonal [11] orientation |24j . 

An analogous effect of orientational pinning has also been seen in exper- 
iments on YBCO superconductors with twin boundaries |45j . In this case, 
due to the correlated nature of the disorder, the direction for easy flow is the 
direction of the twins. A similar effect of horns in the parametric voltage 
curves are therefore observed in the direction corresponding to the twins. 
Also transport measurements when the sample is driven at an angle with 
respect to the twin show a large transverse voltage. 

It is interesting to compare with the angle-dependent transverse voltage 
calculated for d-wave superconductors Also in this case, the transverse 
voltage vanishes only in the [10] and [11] directions. However the 0t vs. cj) 
curves are smooth in this case, since tan#t oc sin4^ |H4j . This is because 
there is no pinning and the transverse voltage is caused only by the intrinsic 
nature of the d-wave ground state. On the other hand, the breaking of 
rotational symmetry studied here is induced by the pinning potential, and 
it results in non-smooth responses like "horned" parametric voltage curves, 
critical angles, transverse critical currents, etc. 

In superconductors with a square array of pinning centers, typically the 
pins are of circular shape and the size of the pins is much smaller than the 
distance between pinning sites [1611181 ITU] . In this case, the pinning barriers 
that vortices find for motion are the same in many directions. Therefore 
it is possible to have orientational pinning in many of the square lattice 
symmetry directions. This explains the rich structure of a Devil's staircase 
observed recently in the simulations of Reichhardt and Nori, where each 
plateau corresponds to orientational pinning in each of the several possible 
directions for orientational pinning. This interesting behavior is not possible 
in JJA, however, since the egg-carton pinning potential corresponds to the 
situation of square-shaped pinning centers with the pin size equal to the 
interpin distance. In this case the only possible directions for orientational 
pinning are the [10] and [01], as we have seen here. 

It is worth noting that many experiments in JJA in the past have been 
done in samples with a diagonal bias. For example, van Wees et al. |38j have 
observed the existence of a transverse voltage in their measurements, which 
was unexplained. From our finding that the [ll]-direction is unstable against 
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changes in the angle of the bias, we conclude that any small deviation in 
the direction of the flow of current, either due to tiny fabrication defects in 
the busbars or to disorder in the critical currents of the junctions [27], may 
explain their observation. Also Chen et aii.|65| have reported a transverse 
angle in measurements in JJA driven in the diagonal direction. In their 
case the effect is antisymmetric against a change in the direction of the 
magnetic field. Since transverse voltages due to the instability of the [11] 
direction are even with the direction the magnetic field, their observation 
can not be explained from our results. This means that they have a Hall 
effect, possibly due to quantum fluctuations. However, they report that their 
transverse voltage had also a component which was even with the magnetic 
field (which was discounted in their computation of the Hall angle). This 
particular spurious contribution can also be attributed to a small deviation 
in the direction of the bias or to disorder effects. From this we conclude that 
in order to study the Hall effect in JJA the most convenient choice would 
be a current bias in the [10] direction where the effect of transverse voltages 
at small deviations in the bias or disorder is minimum. 

When our work was upon completion, new studies of the effect of the 
orientation of the bias in driven square JJA have appeared. Fisher, Stroud 
and Janin |56| have studied some of the effects of the direction of current 
in a fully frustrated JJA (/ = 1/2) at T = 0. In their case a transverse 
critical critical current and the dynamics as a function of I x and I y has 
been described. Their results are in part complementary to our work with 
a single vortex (/ = 1/L 2 ). Yoon, Choi and Kim |57| find differences in 
the IV characteristics of JJA at / = when comparing the parallel current 
bias with the diagonal current bias. Their results are in agreement with our 
results. 

It is worth to mention that very recently an experiment in another sys- 
tems with periodic pinning potentials, in superconductors with a square 
antidots arrays, evidence of guided vortex motion was observed jHS]- They 
showed that the pinning landscape provided by the square antidot lattice 
influences the vortex motion, given place to an anisotropic motion, temper- 
ature dependent, confirming our results of Sec. 1.3.5. 

In summary, in Sec. 1.3.5 of this review we have considered the dynamics 
of a single vortex in a square JJA and diluted vortex lattices. We were able 
to characterize in detail the orientational pinning and breaking of rotational 
symmetry in this case. Furthermore, with the results of the RSJ numerical 
calculation we were able to reproduce and interpret most of experimental 
measurements for a quasi-diagonal bias [21]. Regarding the moving vortex 
lattice different dynamical phases in a JJA, as a function of temperature 
and current, we also show here that characteristics of the breaking of rota- 
tional invariance, orientational pinning and transverse voltages depend on 
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the dynamical phase under consideration as well as on the disorder in the 
Josephson couplings [21] • Therefore we shown that orientational pinning is 
a useful tool to characterize the non equilibrium vortex dynamics regimes 
in JJA, both numerically and experimentally. 

1.4 VORTEX DYNAMICS IN FULLY 
FRUSTRATED JJA 

In this section we show the study of the fully frustrated JJA (/ = 1/2) 
for system sizes of L x L junctions, with L = 8, 16, 24, 32, 48, 64, 128. In 
the absence of external currents, we find an equilibrium phase transition 
at T c = 0.45 which, within a resolution of AT = 0.005, corresponds to a 
simultaneous (or very close) breaking of the U(l) and the Z2 symmetries. 
Here we will analyze the possible occurrence of these transitions as a function 
of temperature when the JJA is driven by currents well above the zero 
temperature critical current I > I c q = (y/2 — l)Io as 0.414Jo- 

1.4.1 U(l) symmetry and transverse superconductivity 

In the driven JJA superconducting coherence can only be defined in the 
direction transverse to the bias current |23| 1251 150 j . We calculate the trans- 
verse helicity modulus as was shown in Sec. 1.2.3. We find that T x is 
finite at low T and vanishes at a temperature Tyn^I). In Fig ll.25f a) we 
show the behavior of T X (T) for a current / = 0.9 in a 64 x 64 JJA. This 
transition is reversible: we obtain the same behavior when decreasing T 
from a random configuration at T = 1 and when increasing T from an or- 
dered state at T = 0, see Fig ll.25f a). Transverse superconductivity can be 
measured when a small current It r is applied perpendicular to the driving 
current: we find a vanishingly small transverse voltage Vt r below Tu^(I), 
as we found before in |23| I25j for / = 1/25. We obtain the voltage in the 
^-direction as the time average = (da^it) / 'dt) (normalized by RnIq); 
longitudinal voltage is V = V y and transverse voltage is Vt r = V x . In 
Fig ll^Sf b) we see that the transverse resistance Vt r /Itr is negligibly small 
for T < Tu(\\ and starts to rise near the transition. The equilibrium U(l) 
transition (at / = 0, I = 0, Kosterlitz-Thouless) is characterized by the 
unbinding of vort ex- ant i vortex pairs above T c . We calculate the density n va 
of vortex-antivortex excitations above checkerboard vortex configuration as 
2n m = (|6(n)|) — /, where the vorticity at the plaquette n (associated to the 
site n) is b(h) = — A M x nint[0 M (n)/27r]. We see in Fig |1.25r c) that n va rises 
near Moreover, the transverse resistivity above is Vtr/Itr oc n va . 

In Fig |1.26l shows T x for sizes L = 32, 48, 64, 128, we see that a transition 
temperature can be defined independently of lattice size. 



40 



V.I. Marconi, D. Dominguez 



. 1 1 1 1 1 1 1 1 1 1 1 1 l" 




• 

i 


in ^ <»* 
* 

T u(i) ( b > 




i i i i i i i i i i i i i _ 

' (cl 



0.00 0.10 0.20 0.30 

T 



Figure 1.25: For large currents I > I c (0): breaking of the £7(1) symmetry. 
Example showed for / = 0.9 and system size 64 x 64. (a) Helicity Modulus 
T x vs temperature T (• increasing T, decreasing T). (b) Transverse 
voltage for a small transverse current, l% r = 0.1, vs T. (c) Vortex-antivortex 
pairs density, n va vs. T. 



1.4.2 Z 2 symmetry 

Since the ground state is a checkerboard pattern of vortices, we define the 
"staggered magnetization" as m s (n, t) = (—l) nx+riy [2b(n x ,n y ,t) — 1] and 
m s (t) = (1/L 2 ) Tn s (h, t). At T = 0, I = there are two degener- 
ate configurations with m s = dbl. Above the T = critical current I c q 
the checkerboard state moves as a rigid structure and m s (t) changes sign 
periodically with time. Therefore we define the chiral order parameter as 
X = {m^(t)). We start the simulation at T = with an ordered checkerboard 
state driven by a current I > I c q and then we increase slowly the tempera- 
ture. We obtain that the chirality parameter vanishes at a temperature Tz 2 , 
which is smaller than as can be seen in Fig ll.27r a) for / = 0.9. This 
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Figure 1.26: Helicity Modulus T x vs temperature T for / = 0.9 and system 
size 64 x 64, increasing T. Size effect for L = 32, 48, 64, 128. 



transition is confirmed by the size analysis shown in the inset of Fig ll.27f a): 
for T < Tz 2 the chirality x is independent of size, while for T > Tz 2 we see 
that x ~ l/L 2 - As it is shown in Fig |1.27f b^. the longitudinal voltage V 
has a sharp increase at Tz 2 , which could be easily detected experimentally. 
The excitations that characterize the Z2 transition are domain walls that 
separate domains with different signs of m s . The length of domain walls in 
the direction \i is given by £ M = (2/L 2 ) ^ fi (6(n)6(n + z/)), with v _L We 
find that for / > and T > the domains are anisotropic, with the domain 
walls being longer in the direction perpendicular to the current {C x > C y ) 
and the domain anisotropy C x /C y increasing with /. In Fig ll.27f c) we show 
the dependence of with temperature for / = 0.9. At T = there are 
no domain walls, since the initial condition is the checkerboard state, and 
the domain wall length grows with T, showing a sharp increase at Tz 2 - The 
domain anisotropy C x /C y is shown in the inset of Fie H^f c): it has a clear 
jump at the transition in Tz 2 while for T 3> Tz 2 the domains tend to be less 
anisotropic. When decreasing temperature from a random configuration at 
T = 1, an important number of domain walls along the x direction remain 
frozen below Tz 2 - C x tends to a finite value when T — > and the domain 
anisotropy tends to diverge when cooling down. This leads to a strong hys- 
teresis in the voltage V at Tz 2 (see Fig |1.27f b)) since the extra domain walls 
increase dissipation |39[ I67j. This low T state with frozen-in domain walls 
is ordered along the x-direction (i.e. perpendicular to I) but is disordered 
along the y direction which gives x ~ 0- We define the Z2 order parameter 
in the x direction as Xx = (0-/L) J2n y i( l / L ) J2n x m s (n x ,n y ,t)} 2 } and x y is 
defined analogously. We see in Fig ll.27l^ a) that, when cooling down from 
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Figure 1.27: For a large currents, / > 7 C (0): breaking of the Z2 symmetry. 
Example shown for a dc current / = 0.9 and system size 128 x 128. Results 
for increasing T (•) and decresing T (<})• (a) Chiral order parameter x vs 
T and one-dimensional order parameter Xx vs. T. Inset: size effect for x 
for L = 8, 16, 32, 48, 64, 128. (b) Longitudinal voltage V vs. T. (c) Domain 
wall lengths C x and C y vs. T. Inset: domain anisotropy jC x /C y vs. T. 



high T, Xx vanishes as Xx ~ V-^ f° r ^ > ^2 has stronger size effects 
than x) an d becomes finite for T < Tz 2 , while Xy ~ for any T. Therefore, 
depending on the history, there are two kinds of high current steady states 
with broken Z2 symmetry at low T. One state has mostly the checkerboard 
structure (x 7^ 0) with few very anisotropic domains. It can be obtained ex- 
perimentally by cooling down at zero drive and then increasing /. The other 
steady state is ordered in the direction perpendicular to / (x x 7^ 0, Xy = 0) 
with several domain walls along the x direction. These domain walls move 
in the direction parallel to / (via the motion of vortices perpendicular to /) 
giving an additional dissipation. This state can be obtained experimentally 
by cooling down with a fixed /. 

The two steady states have also different critical currents as can be 
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Figure 1.28: Current-voltage hysteresis for T = 0.02 shown for (a) voltage 
vs. / and (b) chiral order parameter vs. /. Increasing current from the 
checkerboard state (•) and decreasing current from a random state at large 
I > Iz 2 (CO- Insets: snapshots of the staggered magnetization m s (n, t): 
almost perfect ordered state obtained increasing current (right), and ordered 
state with £ x -domain walls obtained decreasing current (left). 



observed in the low T current-voltage (IV) characteristics. In Fig ll^Sf a) we 
show the IV curve for T = 0.02 and in Fig ll^Hf b) the corresponding \ vs - I 
curve. When increasing / from the 1 = equilibrium state, we find a critical 
current I C 2 (T) , which in the limit of T = tends to I c o = y/2 — 1 as found 
analitically and in simulations with PBC |36| I68| I69j . Near I C 2 the order 
parameter x has a minimum and rapidly increases with /. The driven state 
is an ordered state similar to the one shown in Fig ll.27f a) (see right inset). 
At a higher current Iz 2 there is a sharp drop of x which corresponds to the 
crossing of the Tz 2 (I) line (see Fig ll.29)) and the Z2 order is lost. If we now 
decrease the current either from the disordered state at / > Iz 2 or from a 
random initial configuration or from a configuration cooled down at a fixed 
I > J C 2, we obtain the steady state with domain walls along the x direction 
and x ~ 0) Xx 7^ (see left inset). This state has a higher voltage and pins at 
a lower critical current I c i(T), which has the T = limit I c i(T — > 0) = 0.35. 
It has been shown that open boundary conditions can nucleate domain 
walls leading to the critical current l c i(0) = 0.35 usually found in open 
boundary T = simulations |39l 1461 1471 I67j . Also a moving state with 
parallel domain walls (as in the insets of Fig ll.28f a)) has been found by 
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Gr0nbech- Jensen et al. 70 in / = 1/2 JJA with open boundaries and 
Marino and Halsey [JTj have shown that the high current states of frustrated 
JJA can have moving domain walls. We have studied the effect of open 
boundaries in the direction of I, in the direction perpendicular to I and in 
both directions. They differ mainly in the T = critical current and IV 
curve, for finite T there are small differences in the detailed shape of the 
hysteresis in critical current. In all the cases the two high current steady 
states are observed at finite T with the same history dependence. Also, we 
find that the density of frozen C x domain walls depends on cooling rate and 
decreases with system size. 

1.4.3 Summary 

In summary, we have obtained the current-temperature phase diagram of 
the driven fully frustated XY model, which is shown in Fig ll. 29f a) for low 
currents and in Fig |1.29f b) for high currents. At high currents the breaking 
of the U(l) and the Z2 symmetries occurs at well separated temperatures, 
with Tz 2 < Tjjm . The low temperature regime T < Tz 2 (I) has bistability 
with two possible steady states (see insets in Fig ll .29|) and history dependent 
IV curves. The different transitions could be observed experimentally with 
measurements of the transverse and longitudinal voltage. 

1.5 CONCLUSIONS 

We presented a detailed review of our previous numerical studies on non 
equilibrium vortex dynamics in JJA driven by a dc current. Dynamical 
phase diagrams for different magnetic field corresponding to / = 1/25 and 
/ = 1/2, current directions (orientational pinning) and varying temperature 
were obtained and shown in Fig ll. 24l and Fig ll. 29l 

In conclusion for diluted vortex lattices we find that for low currents 
the depinning transition of the VL and the later melting of the moving VL 
become separated with Tm{I) > T P {I). For large currents we find that the 
moving VL has a finite transverse critical current and therefore transverse 
superconducting coherence. In this case it is possible to define a transverse 
depinning transition at a temperature Tf r , and a later melting transition 
of the moving VL at Tm(I). This transverse depinning transition could be 
easily studied in controlled experiments in Josephson junction arrays, both 
with transport measurements and with inductive coil measurements of the 
transverse helicity modulus. This proposed experiment were later performed 
|61j and the anisotropy of moving vortex phases was confirmed. 

Regarding fully frustrated JJA we obtain the current-temperature phase 
diagram. At high currents (Fig ll.297 b)) we find Tz 2 < ^c/(i)j which is in 
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Figure 1.29: Current-temperature phase diagrams: (a) for low currents, I < 
2c(0) and (b) for high currents, I > I c (0). I c i(T) (0) and I C 2(T) (•) lines are 
obtained from hysteresis in IV curves as well as from hysteresis in T x (T) and 
X {T) curves. T U(1) (I) line obtained from T X (T) and V tr (T) (A). T Z2 (I) line 
obtained from x{T), V(T) and C x /C y (T) (•). The dashed line corresponds 
to the IV curve of previous Fig. (T = 0.02). Insets: Two different steady 
states (staggered magnetization snapshots) found below Tz 2 (I) line, at T = 
0.02 and / = 0.9: (i) warming up, almost a perfect checkerboard pattern, or 
(ii) cooling down, ordered state with /^-domain walls. 
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contrast with the equilibrium result of Tz 2 — Tu(i) [1 3 5]- Interestingly, the 
frustrated XY model with modulated anisotropic couplings also has Tz 2 < 
Tu(i) > meaning that the anisotropy induced by the current may provide 
a similar effect. It is clear that the / = 1/2 case has a strong pinning effect 
(with I c q = 0.414) when compared to the dilute case of / = 1/25 (with I c q = 
0.11) jlHlEnj. In fact, the transverse depinning temperature Tyix) is one 
order of magnitude higher for / = 1/2 with respects to / = 1/25 [231- The 
driving current weakens the effect of pinning and thus Tu(i) increases with /. 
A similar effect gives a Tz 2 growing with I just above I c q. However for larger 
currents (near the Josephson current Iq) Xz 2 (I) starts to decrease with /, 
with the limit Tz 2 (I — ► oo) —* 0. This is because a driving current increases 
the density of domain walls (an effect already mentioned in [SHj) destroying 
the Ising order for / 3> I c o- Moreover, we find that the ordered region 
in T < Tz 2 (I) has bistability with two possible steady states and history 
dependent IV curves. The different transitions here characterized could be 
obtained experimentally by measurements of transverse and longitudinal 
voltage. With our numerical results on highly driven fully frustrated JJA 
we expect to motivate an experimental confirmation of the clear separation 
of the transition temperatures, Tz 2 < ?t/(i)- 



Appendix A: Periodic boundary conditions 

Phases 

We want to obtain the periodic boundary condition (PBC) for supercon- 
ducting phases 6(n x ,n y ). In general we can write the PBC as: 

#(n + L x ) = 0(n)+u x (n) 

#(n + L y ) = d(n)+u y (n) , (1.28) 

where L x = (L x , 0) and L y = (0, L y ). Taking into account that all variables 
should be independent of the order of subsequent global translations and 
that the phases are defined except for an addition of 2ttI, we are lead to the 
consistency condition 

u x (n) + u y (n + L x ) = u^n + L y ) + %(n) + 2nl , (1-29) 

Therefore, in order to specify the periodic boundary conditions we have to 
give the functions u x ,u y , which will depend on the gauge for the vector 
potential A. 

The periodic boundary condition for the phases can be deduced by re- 
quiring that all physical quantities are invariant after a translation in the 
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lattice size. The supercurrents S^n) = sin[A A1 0(n) — yl M (n)] should satisfy: 

5 M (n + L x ) = S^n) 

S^n + Ly) = S^n). (1.30) 

This implies that the gauge invariant phase difference #u(n) = A^8(n) — 
A^(n) should satisfy 

^(n + L x ) = ^(n) + 2vr/ 

6» M (n + L y ) = ^(n)+27rZ', (1.31) 

with any integer. This condition leads to 

A^(n) = 2tvI + A^in + L x ) — A^in) 

A M% (n) = 2vrr + ^(n + L y ) - A M (n) . (1.32) 

We can choose the solution with I = I' = 0. In the Landau gauge, A x is a 
linear function of n y and A y is a linear function of n x . Taking the origin 
such that A^n = 0) = 0, we obtain for the Landau gauge: 

tta.(n) = u x (0) + A y (L x )n y 

u y (n) = u y (0) + A x (L y )n x . (1.33) 

The consistency condition ()1.29j) requires 

A y (L x )Ly - A x (L y )L x = 2irl . (1.34) 

The term in the left side is equal to the total flux 2nfL x Ly, therefore (|1.34|) 
is equivalent to flux quantization, giving fL x L y = N v , with N v = I the 
number of vortices. 

If we take u x (0) = u y (0) = 0, we obtain for the PBC 

8(n x + L x ,n y ) = 0(n x , n y ) + A y (Jj x )n y 

6(n x ,n y + L y ) = 9(n x ,n y ) + A x (L y )n x . (1.35) 

A particular choice can be the gauge with A x (n) = —2-Kfn y , A y (n) = 0, 
which leads to Eq. 11.81 

External currents and electric fields 

In the presence of external currents or voltages the periodic boundary condi- 
tions have to be reconsidered. In this case it is possible to have £ E • dl ^ 
in a path that encloses all the sample either in the x or the y direction. 
Therefore, in a closed path we have to consider the Faraday's law 



/ E . d i=-i* 

J c dt 
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The two-dimensional sample with PBC can be thought as the surface of 
a torus in three-dimensions. The closed paths we are considering are the 
two pathes that encircle the torus. The electric field is not a gradient of a 
potential, it is now given by E = — W — One possible solution is to 

consider a vector potential 



A(r,t) = A (r) - 5(t), 

for which 



H cxt = VxA 
Ida 

Etotal = - -77 

C at 

In our case, we take the adimensional vector potential as: 

A li (n,t) = A° tl (n)-a li (t), (1.36) 

with A°(n) in the Landau gauge (A®(n) = —2irfn y , Ay(n) = 0). Therefore 
the gauge invariant phase is: 

0„(n, t) = A M #(n, t) - Aj(n) + <*„(*). 

Then a^ acts as a global time-dependent phase in the \i direction. 

In the normalized units used in this paper, the electric field in the link 
defined by the junction n, fi is 

E „(„) = -A„V(n) - 
where the electrostatic potential is V(n) = — de ^ . Therefore we have 



(it 

d6(n) da 



^(n) = A M -^ + 



fi 



dt dt 

The average electric field in the /x direction is: 

da 



where we have used the fact that ^ n A^ de ^ = (which is the discrete 
equivalent of §W • dl = 0). The current in the link n, \x is, in normalized 
units,: 

/ M (n)=£ M (n)+^(n), 
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with S^(n) = S^(n) + i]^(n,t). Therefore the average current in the fi 
direction is 



T 



av 

^ L X Ly 



There are two cases to consider: (i) external current source: the external 
current is given and the total voltage fluctuates and (ii) external voltage 
source: the external voltage is given and the total current fluctuates. 

(i) External current source: the average current in fi direction is fixed by 



7"ext 

quantity given by: 



the external current: 1^ = 1° . The average electric field is a fluctuating 



In this case, the ctu(t) is a dynamical variable, and its time evolution is given 
byEq.[071 

(ii) External voltage source: the average electric field in the [i direction 
is fixed by the external electric field: E^ = E^ 1 . Therefore, now the is 
given by the external source: 



and OLfiii) = Ef^t, if E c * 1 is time-independent. The average current is now 
a fluctuating quantity given by: 



/- = jB ext + _2_^ (n) . 



ii 



Let us see how the PBC are affected by a change of gauge. The gauge 
transformations are the following: 

0(n) -> 0'(n) = 0(n) + 0(n) 
^(n)^4(n) = A M (n) + A^(n) 

(n) 



V(n) -> y'(n) = ^(n)- ^ 



An interesting choice is: 



0(n) = a • n = i 
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In this gauge we have ^(n) = A„0(n) — -A^(n) and the PBC for the phases 



is: 



0(n + L x ) = 6(n) + A y (Lx)n y + a x L x 

6(n + L y ) = e(n) + A x (L y )n x + a y L y , (1.38) 

and for the voltages: 

ft CM 

V(n + L„) = V(n) - -^L, = V(n) - ^ (1.39) 

The equations of motion in this gauge are: 
d9(n) ^ ^ a 



dt ' dt 

H n' 

The periodic boundary conditions with a fixed external current using 
Eq. 11.371 was used previously in Ref. [72] for / = (it was called a "fluc- 
tuating twist boundary condition") and in Ref. |35j for f ^ 0. Also, the 
periodic boundary conditions in the gauge of Eqs. 11.381 and 11.391 were used 
in Ref. [HI] for a time dependent d-wave Ginzburg-Landau model. 



Helicity modulus 

The helicitity modulus expresses the "rigidity" of the system with respect 
to an applied "twist" in the periodic boundary conditions. The twist is 
defined as a phase change of L^k^ between the two opposite boundaries 
which are connected through the PBC in the fi direction 

9(n + L^=e(n)+L^. 

The helicity modulus is obtained from the free energy F(T, k^) as: 

1 d 2 F(T,k^) 
^ L 2 dk 2 

It is clear from Eg. 1 1.381 that fc„ = a». Then, in order to evaluate the helicity 
modulus, must be set to zero. This means that the helicity modulus can 
not be calculated in the direction in which there is an applied current, since 
it gives a fluctuating twist k^t). 



(1.41) 
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Appendix B: Algorithm 

The core of the numerical calculation is to invert the Eq. 11.101 This means 
to solve a discrete Poisson equation of the form, 

Aj/(n) = d(n) (1.42) 

with the periodic boundary conditions 

/(n + L x ) = /(n) 

/(n + L y ) = /(n). (1.43) 

The linear system of L x L y equations of (jl.42j) is singular. Physically, the 
reason is that in Ea. ([1.10|) . /(n) = d9(n)/dt corresponds to a voltage, which 
is defined except for a constant. We choose the voltage reference such that 
it has zero mean: 

E^ n ) = ' (1-44) 

n 

other choices are also possible (like, for example, fixing /(no) = at a given 
site no). 

The method we use to invert the Eq. (|1.42|) is based on the Fourier 
Accelerated and Cyclic Reduction (FACR) algorithm 73 J. In this case, we 
take first the discrete Fourier transform in the x direction: 

d(k x , n y ) = } j d(n x , n y )e L * , (1-45) 

Hi 

which with a fast Fourier algorithm takes a computation time of order 
L x \ogL x . This leads to the following equation: 

£k x f(k x ,n y ) - f(k x ,n y - 1) - f(k x ,n y + 1) = d(k x ,n y ) , (1.46) 

with ek x = 4 — 2cos^^ and with boundary condition f(k x ,L y + 1) = 

f(k x ,l). This is a cyclic tridiagonal equation which can be solved with a 
simple LU decomposition algorithm in a computation time of order L y |73j . 
In this way the f(k x , n y ) is obtained from Q1.46JI . Finally we take the inverse 
Fourier transform to obtain 

f(n x , ny) = — f(k x ,n y )e . (1.47) 

X h 

This algorithm takes a computation time which is of order L a; L 3/ (logL :r + ^4) 
with the constant A ~ 1. This is faster than the two dimensional Fourier 
transform method used by Eikmans and van Himbergen |48| . which takes a 
computation time of order L x L y (logL x + \ogL y ). 
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